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Abstract 

We present an extension of the Allen-Cahn/Cahn-Hilliard system which incorpo- 
rates a geometrically linear ansatz for the elastic energy of the precipitates. The model 
contains both the elastic AUen-Cahn system and the elastic Cahn-Hilliard system as 
special cases and accounts for the microstructures on the microscopic scale. We prove 
the existence of weak solutions to the new model for a general class of energy func- 
tionals. We then give several examples of functionals that belong to this class. This 
includes the energy of geometrically linear elastic materials for D < 3. Moreover we 
show this for D = 3 in the setting of scalar-valued deformations, which corresponds to 
the case of anti-plane shear. All this is based on explicit formulas for relaxed energy 
functionals newly derived in this article for D = 1 and D = 3. In these cases we can 
also prove uniqueness of the weak solutions. 

1 Introduction 

In this article we study the Allen-Cahn/Cahn-Hilliard model (AC-CH model for short), 
which combines and extends two famous diffuse interface models: the AUen-Cahn equation 
and the Cahn-Hilliard equation. Both have been applied successfully to model segregation, 
precipitation and phase change phenomena in alloys and liquid mixtures in materials science, 
geology, physics, and biology, among others. 

The AC-CH system was first introduced in [Tl] , and its mathematical properties have been 
studied extensively, see, e.g., [31 [37] and references therein. We study here an extension of 
this model with a particular ansatz for the elastic energy. 

The AUen-Cahn equation was first introduced without elasticity in [T]; and the Cahn-Hilliard 
equation was first introduced without elasticity in |12j . The AUen-Cahn equation is a 
second order partial differential equation of Ginzburg-Landau type for an unconserved order- 
parameter and thus can be used to model segregation and precipitation in solids, or other 
more general situations where a reordering of the crystal lattice occurs. Conversely, the 
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Cahn-Hilliard equation is a fourth-order partial differential equation for a conserved order 
parameter striving to model phase change phenomena where one physical quantity like the 
volume of the phases is preserved. 

The Allen-Cahn system with linear elasticity was studied before in [3], the Cahn-Hilliard 
system with linear elasticity in [3T] , [53] and [13] . An extension of the Cahn-Hilliard system 
with geometrically linear elasticity valid for single crystals was recently found in |7] for 
D<2. 

Except for the latter work, elastic effects due to small scale microstructures within the 
phases have been neglected. Here we treat the combined AC-CH model with elasticity in 
D < 3 dimensions. We state the new extended model in ([J)-© in Section |21 This model 
provides a basis for the generalization of further existing isothermal diffuse interface models. 
Our definition of the extended model does not require any special assumptions (except 
regarding the regularity) on the stored energy functional W. However, for the proof of 
existence of weak solutions, we require the Assumption (A) phrased in Section [3] below. 
We show that Assumption (A) is satisfied in the following cases: 

(i) Materials which follow the linear theory of elasticity developed by Eshelby [T^ in the 
context of elastic inclusions and inhomogencities (see a remark after Assumption (A)); 

(ii) materials in < 2 which are well-described by a geometrically linear theory of elas- 
ticity (Theorem 13 for £> = 2, cf. also [TU]): 

(iii) materials in D = 3 which are well-described by a geometrically linear theory of elas- 
ticity, where the deformations are assumed to be scalar-valued functions. This corre- 
sponds to the situation of anti-plane shear, see below for details and cf. Theorem |5l 

The geometrically linear theory of elasticity allows for fine microstructures within the phases 
of the elastic materials taken into account. This becomes apparent in the explicit formulas 
for the relaxed energy functionals, which are, together with formulas for their derivatives, the 
basis of our proof of Theorem [5] In Section |4| we recall an explicit formula in D = 2 derived 
by Chenchiah and Bhattacharya [15]. Moreover we prove explicit formulas for D = 1, and 
for D = 3 in the setting of scalar- valued deformations. (In [TS] partial results for = 3 in 
the vectorial, i.e. the non-scalar, setting are shown.) 

The coupling to elasticity changes significantly the morphology of the precipitates and the 
coarsening patterns, see, e.g., the classification in [5D]. For an intuitive picture of the 
microstructures taken into account we refer the reader to Figures [1] and [2] Figure [1] 
displays exemplary macroscopic phases of a segregated material with phases d being almost 
one (dark gray) and almost zero (light gray). In Figure [2] we show the corresponding 
microscopic length scale featuring microstructures; the left is a blow-up of a small region 
within the phase with d being almost one and the right is a blow-up of a very small region 
within the phase with d being almost zero. On this microscopic scale where we assume that 
it is sufficient to treat the elastic energies of single crystals, fine microstructures in the form 
of laminates occur, see, e.g., [301 [^ . 

The methods developed here apply generally to any established phase change and segre- 
gation model provided the temperature is conserved. (For non-isothermal settings, the 
validity of the second law of thermodynamics requires additional corrections which are not 
studied here.) For further discussions of our model and the analytical results we refer to 
the conclusions in Section [6] 
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Figure 1: Macroscopic phases of a segregated material with phases d being almost one (dark 
gray) and almost zero (light gray) . In Figure [2] blow-ups of the circular regions are shown. 



2 The AC-CH model and extensions 

Throughout this paper, let fl C M.^ for D > 1 be a bounded domain with Lipschitz boundary 
which serves as an (unstressed) reference configuration. For a stop time T > 0, let 51 7- :~ 
il X (0,T) denote the space-time cylinder. To the AUen-Cahn/Cahn-Hilliard system, first 
derived in |14) . we add elasticity, possibly respecting the lamination microstructure, by 
introducing the system 

dta = Adiv(Af(a,6)V^), (1) 

dtb = -M{a,b)^, (2) 
= div(aeW(a + 6,e(u))), (3) 

where the function a : ilj- — !■ is a conserved order parameter, typically a concentration, 
b : Qj- —5- M.Q is an unconserved order-parameter, specifying the reordering of the underlying 
lattice, M{a, b) denotes the positive semi-definite mobility tensor, A > is a small constant 
determining the interfacial thickness. The choice of W determines whether lamination 
microstructure occurs in the system, see (|2.13|) . 

By u : O — )- we describe the displacement field, such that a material point x in the 
undeformed body is at x' = x + u(x) after the deformation. Then the (linearized) strain 
tensor is defined by 

e(u) :=i(Vu + Vu*), (4) 

where A* denotes the transpose of a matrix A £ R-^^^. As usual, • stands for the inner 
product in R^, that is u- v = X^^Li '^i^i: for A, B £ R^^^ we denote the inner product 
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Figure 2: Microstructurcs with d being almost one (left) and almost zero (right). The black 
part corresponds to d = 1, while d = in the white regions. The average of d equals d. 



pDxD 



by 



A:B := tr(A*B) 



D 



Moreover, \A\ := \JA:A for A € R"^^-^ is the Frobenius norm. We denote the symmetric 
matrices in M^^-° by R^y^J^ . 

The functional 'W{a + b,e{u)) represents the stored elastic energy density. We will choose it 
either according to the linear ansatz by Eshelby, W = Wun, or as a single crystal composite 
lamination energy, W = W. The definition of W is given in the coming section \2A\ 
The linear theory by Eshelby [T9| developed in the context of elastic inclusions and inho- 
mogeneities, can be summarized in the following ansatz for the elastic energy 
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WUd, e) :=-(£- sid)) ■ C{d){e - e{d)) 



(5) 



for all £ e M.^yll^, d-~a 



de with a constant e € R£^^. 



Dy.D ^ „ , 5^ and £(d) 
By C{d) we denote the symmetric, positive definite and concentration-dependent elasticity 
tensor of the system that maps symmetric tensors in MP^^ to themselves. 
The system ©-(jSj is completed with the definition of the free energy 



F(a, &, u) := j V(a, b) + - {\Va\^ + \Vh\'') + W{a + 6, £(u)) + W.^Mn)) Ax, (6) 
n 

see [H], where ■0(0, h) is the free energy density assumed to be 

— yg{a + h) + g(a - b)j + Kia(l — a) — K^b^ , 
g{s) s In s + (1 — s) ln(l — s) 



ij{a,b) 



(7) 



for scalars ki, K2 > 0. The term ^[^(a + b) + g{a — b)] in ([7]) defines the entropic part of 
the free energy, given in the canonical Bernoulli form for perfect mixing, and > is the 
constant temperature. 
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The functional Woxt(£) in © represents energy effects due to applied forces. In the absence 
of body forces, the work necessary to transform the undeformed body f2 into a state with 
displacement u is then 

- / u • CToxtndS' = - / Vu : cToxt dx = - e(u) : (Text da;, 
Jdn J J 

Q n 

where we use that the applied stress (Text is constant and symmetric. Consequently, 

W^cxt(e) = -e : acxt (8) 

is the energy density of the applied outer forces. The system ([T])-([3]) has to be solved in flj- 
subjcct to the initial conditions 

a(< = 0) = ao, 6(t = 0) = 60 in 

for given functions oq, 60 ^ ^ IK. subject to the Neumann boundary conditions for a, the 
no-flux boundary conditions, and the equilibrium condition for applied forces 

Va • n = 0, J(a, b,u) = 0, a ■ n = (Text ■ n on dil, t > 0. (9) 

Here, a := d^yV{a + b, e{u)) defines the stress. 

In ([5]), n is the unit outer normal to dfl. For simplicity, body forces are neglected and it is 
assumed that the boundary tractions are dead loads given by a constant symmetric tensor 
foxt- By J we denote the mass flux, given by 

dF 

J(a,6,u) := -M(a,bWiJ = -M(a,bW — (a,b,u), 

oa 

with ^ :~ ^ the chemical potential. 

The valid parameter range of a and b is, see Theorem [TJ 

0<a + 5<l, 0<a-&<l. (10) 

The inequalities are strict unless (a, b) = (0, 0) or (a, b) = (1, 0). The system ©-(IS]) includes 
as special case the elastic Cahn-Hilliard system (setting 6=0, |21) ) 

/ f)F\ 
dta = Adivf M(a)V—, 
\ oa / 



with 



div(5eW(a,e(u))) 



F{a,u) := J 7/.(a) + ^|Vap + W(a,£(u)) + Wext(£(u))da;, 

n 

■d 

il>{a) := — alna + (1 — a) ln(l — a) + Kia(l — a) 
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and the boundary and initial conditions correspondingly to above. Moreover the system 
(IH)^® includes as special case the elastic Allen-Cahn equations (setting a = ^, [9]) which 
for rcscaled h with < & < 1 read 

dF 



= div(9e>V(6,e(u))), 



with 



F{h,n) := j V(6) + ^|V6|2 + W(6,£(u)) + iyext(e(u))da:, 



m ■.= -{9{b)+g{l-b))+^~n,b^ 



and g as well as the boundary and initial conditions are as above. 

The system (HI)-® is exemplary for an isothermal model that exhibits simultaneous ordering 
and phase transitions. Equation ([T|) is a diffusion law for a governed by the flux J and states 
the conservation of mass in fl. Equation ([2]) is a simple gradient flow in the descent direction 
— ^ . Equation ([3]) is a consequence of Newton's second law under the additional assumption 
that the acceleration duu originally appearing on the left hand side can be neglected (this 
can be proved formally by a scaling argument and formally matched asymptotics). The 
vector equation ® serves to determine the unknown displacement u. 

Remark 1. The equations ([2P~([3i can he generalized to vector-valued mappings a, b. This 
allows to study situations with more than two phases present. To fix ideas and for the sake 
of a clear presentation, we restrict ourselves throughout this paper to scalars a and b. 

Remark 2. The equations ([2P-(0i with boundary conditions and and ([7]) comply 
with the second law of thermodynamics, which in case of isothermal conditions reads for a 
closed system 

dtF{a{t),b{t),u{t)) < 0. 
This inequality can be verified by direct inspection similar to the calculations in J^. 

In the proof of Theorem[5]we apply the following explicit formulation of ([I])"® with constant 
mobihty M = 1 



dta = AA 



1? dW 

-{g'{a + b)+g'{a-b)) + Ki(l-2a) + —{a + b,e{u)) - Aa 



(2.1') 



I? dW 
dtb = XAb +- [g'ia -b)-g\a + b)] + 2K2b-^ia + b,eiu)), (2.2') 

= dw{deWia + b,e{u))). (2.3') 



The subsequent section is devoted to the geometrically linear elastic energy for single crys- 
tals. This is a prerequisite to the discussion of existence and uniqueness results for the 
extensions of the Allen-Cahn/Cahn-Hilliard (AC-CH) model studied in Section [31 
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2.1 The geometrically linear theory of elasticity in single crystals 

Our main objective in this subsection is to study a geometrically linear theory of elasticity 
in the context of isothermal phase transitions. For systematic reasons, we first recall the 
linear ansatz dating back to Eshelby, [12] . As a byproduct of the existence theory proved in 
Section [3] we obtain a new existence result for the AC-CH equations with linear elasticity. 
Subsequently, we introduce the geometrically linear elasticity theory that takes the laminates 
of the material into account. 

In the following wc assume that two phases arc present in the considered material which 
may form microstructures as displayed, e.g., in Figures |3] and |4l We refer to the energy Wi, 
i = 1, 2 of each of the phases as microscopic energy, cf. (|2.12p . and to the energy W{d, e(u)) 
in (|2.13p . which reflects the effective behavior of the system with microstructures, as the 
macroscopic energy. 

To determine the energy W{d, e{u.)) in the geometrically linear theory we need to solve a 
local minimization problem, Eqn. (|2.13p below, which we shall outline now. 
We assume that the volumes occupied by each of the two phases in are measurable sets. 
In particular, if di = d, d2 = 1 ~ d characterize the two phases on the microscale, we have 
di e BV{n; {0, 1}) and c?i +^2 = 1 a.e. in ft. The symbol BV denotes the space of functions 
of bounded variation, see, e.g., [11I3T]. By 

(g) j'J{x)dx — 1^ y gix)dx (2.11) 
o 

wc denote the average of a function g in fl, where \E\ is the ZJ-dimcnsional Lcbcsgue measure 
of a set E. 

Let ef € R^,^^, i = 1, 2, be the stress-free strain (or eigenstrain) of the i-th phase relative 
to the chosen reference configuration and ai be its positive definite elasticity tensor. Then 
the clastic energy density of phase i subject to a strain s is given by 

W,{e) (e-ef) : {e - ef) + (2.12) 

for Wi > 0. In (|2.12p we assumed for simplicity that ai, ef and Wi are constants, independent 
of a: e f2 and the order parameter d. 

Under the assumption that the elastic energy adapts infinitely fast and that the surface 
energy between laminates of the microstructure can be neglected, the effective elastic energy 

is, [n], 

W7(d,e)(y):= inf inf j W{d,e{\i))Ax, d&[Q,l], (2.13) 

de{0,l} 

where we used 

W{d,e) := dWi{e) + (1 - d)W2{e), de {0, 1}. (2.14) 

The definition (|2.13p requires further clarification. Firstly, e = e{u) := ^(Vu + Vu*), and 
instead of integrating over D,, one may integrate over Br{y), the open ball of radius r around 
y S ri, where the mean (•) is now taken over Br{y). By homotopy arguments or by results 
in [18], any r > yields the same value of W{d,e){y), as long as Br{y) C 51. Taking 
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the union of such balls then leads to (|2.13p . Secondly, the infimum over d is the result of 
homogenization subject to the constraint that the volume fraction of the selected phase is 
preset by d{y), see [T71 Chapter 10]. This infimum is taken over functions d £ BV{^1; {0, 1}) 
as explained above ensuring W > 0. This is why (|2.13p is only meaningful for d G [0, 1]. 
The second infimum is taken over functions u G H^{n; M.^) where the condition u.\dn = 
e{y)x has to be read as u(a;) = £{y)x for a.e. x G dVt. This originates from the requirement 
that the functional W thus defined must be quasi-convex, see [TH], and is the result of 
relaxation theory, [TH], [13]: as follows. If for prescribed d = a + h the microscopic elastic 
energy density is denoted by Wdie)^ then 

Wd{e):= inf -fwd{e{n))dx (2.15) 

is the elastic energy density of the material with macroscopic strain e after microstructure 
has formed. For D ~ 2, explicit analytic formulas for W are known, see [15] as well as 
formulas of the partial derivatives of W. These will be recalled in Section |3J 



3 Existence and uniqueness results for the AC-CH sys- 
tem 

The existence of solutions to the AUen-Cahn/Cahn-Hilliard equation without elasticity was 
studied in [11] with the help of a semigroup calculus. Existence and uniqueness of weak 
solutions to the Cahn-Hilliard equation with linear elasticity is proved in |21| , with geomet- 
rically linear elasticity in [7] . Existence and uniqueness of weak solutions to the AUcn-Cahn 
equation with linear elasticity is shown in [3] . 

Subsequently we provide existence and uniqueness results for ([T])-®, where W is an elastic 
energy density satisfying the following assumption {A). 

(A) The elastic energy density W G C^(M x K^^^; M) satisfies the conditions 
(v41) i9eW(c?, •) is strongly monotone uniformly in d, i.e., there exists a constant 
ci > such that for all £i, £2 € K^JJ^^ and all d G M 

id,W{d,S2) - 9eW(d,£l)) : (£2 - ^l) > Ci|£2 - £lp. 

{A2) There exists a constant Ci > such that for all d G M and all e G K^J^^^ 

\W{d,e)\ < Ci(|d|2 + |£|2 + i), 
|W(d,£)| < Ci(M|2 + |£|2 + i), 
\deWid,e)\ < Ci(|d| + |£| + l). 

All constants in this article, unless explicitly stated otherwise, may depend on the material 
parameters ai, 0:2, ej and £^, but are independent of d and e. Condition {A2) states that 
W is convex in e. The problem becomes non-convex through the dependence on d. One 
prominent example satisfying assumption (A) is the clastic energy Wiin in ([5]). In Section [5| 
we prove that also the relaxed energy functional W defined in (|2.13|) satisfies the assumption 
(A). 
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Theorem 1 (Existence of weak solutions). Let the mobility tensor M be positive definite, 
W fulfill (A), ip be given by ^ and the initial data (ao,&o) satisfy U(^) . 
Then, there exists a solution (a, b, u) to ^\)-^ that satisfies 

fi)a,beC'^'H[0,nLH^)), 

(ii) dta, dtbeL^inr), 

(lii) u g (o,r; H^{n- R^)), 

(iv) The feasible parameter range of {a,b) is given by hlO\) . 

Proof: The statements of the theorem can be proved with the methods developed in for 
an Allen-Cahn system with linear elasticity. We sketch the main steps. 
First we introduce the operator M associated to w i-J- —MAw as a mapping from H^{il) 
to its dual by 

M{w)r]:= J M\/w -yrjdx. (3.16) 

n 

From the Poincarc inequality and the Lax-Milgram theorem (which can be applied since M 
is assumed to be positive definite) we know that M is invertible and we denote its inverse 
by Q, the Green's function. We have 

(MVg/, V77)i2 = (,y,/) for aU e H\n), f e (H^n))' . 

For /i, /2 e {H^{^l)y, we define the inner product 

(/i,/2)Af := (MVe/i,Vg/2)L2 

with the corresponding norm 

II/IU/ yC/TW for/e (i/i(r!))'. 

For a small discrete step size h > 0, chosen such that Th^^ e N, for time steps to € N with 
< m < Th^^, and given values a™^^, 6™^^ £ R, we introduce the discrete free energy 
functional 

F"^^\a, b, u) := F{a, b, u) + ^l|a - a™-'llL + ^11^ - ^""'lll^, (3-17) 

where (in case of to = 1) it holds a'^ = ao, b'^ = bo, the initial values of a and 6. By the direct 
method in the calculus of variations and Assumption (A), it is possible to show that for h 
sufficiently smah, i^™-'' possesses a minimizer (a"", u™) e (n) x (n) x {n; R^). 
This minimizer solves the fully implicit time discretization of ([Ij-©. Next the discrete 
solution is extended affine linearly to (a, 6, u) by setting for t = {rm + (1 — t) (m — l))ft, with 
suitable re [0,1] 

(a,6,u)(t) := T(a",fe",u") + (1 - T)(a™-\ u"^"!). 

The validity of the second law of thermodynamics (cf. Remark [5]) implies that F is non- 
increasing in time. This allows to derive uniform estimates for (a, 6, u). Compactness 
arguments then allow to pass to the limit h\0 and the limit solves (Il])--(l3|). □ 

In general, the uniqueness of solutions to ([I])-® is open. However, we prove it in a special 
case for the linear elastic energy density W = Wun- 
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Theorem 2 (Uniqueness of solutions for linear elasticity). Let W=Wiin be given by the 
material be homogeneous, i.e. the elasticity tensor C be independent of d, and let Af = 1. 
Then the solution (a, 6, u) of Theorem]^ is unique in the spaces stated there. 

Proof: The proof is very similar to the proof of uniqueness for the Cahn-Hilliard system 
but we repeat it here because we later need to modify it. 

Fix to e (0,T). Let {a'' ,b'' ,u''), fc = 1,2 be two pairs of solutions to (|2T|) -(|23 ^ and ©. 
The differences a := — a^, b := b^ — b^, u := — with corresponding difference of the 
chemical potentials ^ :— — :— ^{a^, b^) — §f (a^, b^) solve the weak equations 



J [~adt^ + XV ■ VC] da; dt = 0, 

J [dtbr] + AV6 -Vrj-e-.C (e(u) - (a + b)e) tj] dx dt 



(3.18) 



Or 



Or 



{g'{a' + b') - g'{a' + b') + g'{a^-b^) ~ 9'{a'~b')) V + 2«2&?; 



J C (£(u) - e{a + b)) : e{u) dx di = 



dxdt, (3.19) 
(3.20) 



for every € i^(0,r; H^i^)) n £°°(17r) with dt^, dtrj € L^i^r), i{T) ^ 0, where in 
order to get p.20p we plugged in (u^ — u^)X{^Q to) ^ test function and integrated by parts. 
As a test function in p.lSp wc pick 

^^'^■"1o, ift>to. 



This shows 



J an + AV(^a) • V{dtQa) dx dt = 0. 



(3.21) 



The difference of the chemical potentials fulfils, with the help of 



Kdxdt^ I \-[g\a'+b')^g\a' + b')+g'{a'-b')-g'{a'~b'))(: 

dx dt. 



- 2KiaC + AVa • VC - £ : C(e(u) - (a + b)e)C 
We pick C (a^ - a'^)Xi^o,to)- With (jX^ we obtain 
^lla(io)|iL+ J A|Vap -ae : C(e(u) -£(a + &))da:dt < 



2Kia^ + - 



\g'{a^+b^)-g'{a^ + b^)\ + \g'{a'^-b^)-g'{a'^-b'^)\ \a\dxdt. (3.22) 
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In p.l9p we choose 77 := {b^ — 6^)A'(o,to) as a test function and add the resulting equation 
to ((3^ and use ((3^ . We end up with 

^\H-to)\\M + l\\Hto)\\L-^ + J [A(|Vap + |V&p) +W(a + 5,£(u))] dxdt 
< J 2(Ki|ap + Kalends dt 

\g'ia^ + b^)-g'{a^ + b'')\ + Ig'ia^ -b^) - g'ia^ -b'')\\{\a\ + \b\) dxdt. 



2 

From Thcorcm[l]wc know that the terms g'(a' ± 6*), « = 1, 2 are finite, and g' is Lipschitz 
continuous. Applying first Young's inequality, then Gronwall's inequality, as to € (0, T) was 
arbitrary, we find a = 6 = in Q-j-. This finally yields 



J e(u) : Ce{u)dxdt = 0. 



With Korn's inequality this proves u = in ilj-. 



□ 



4 Explicit formulas for W 

In many situations like the numerical implementation of the extended models, the above 
definition (|2.13p of W is not practical since it is indirect and based on a local minimization. 
For these applications and for direct later use, we collect here some explicit formulas of the 
relaxed energy W for D < 2 and for the scalar case in Z? = 3. 



4.1 The case D = 2 
As shown in [15], it holds 

W{d,e) = diWiiel) + d2W2{e;) + I3*did2 dct(e* - el), (4.23) 

where /?*, el and £2 ^tre defined below. First we need to fix further notations following [l^. 
Let 7* > be given by 

7* := min{7i,72}, (4.24) 

where 7^ is the reciprocal of the largest eigenvalue of ' Ta^ , cki is the elastic modulus 
of laminate i, and the operator T : M.'^y^ — > I^syin given by 

Te = £ - tr(e)Id. 

In [7] a recipe is given for the practical computation of 7*. Here, we only remark that if the 
space groups of the two existing laminates are cubic, it holds 

7* = min{Ci,ii — C1.12, 6*2,11 — 6*2,12, 2Ci,44, 26*2,44}- 
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The first subscript of C denotes here the phase, the other two indices are the coefficients of 
the reduced elasticity tensor in Voigt notation, [28] . 

As shown in [IS], the scalar (3* S [0, 7*] determines the amount of translation of the laminates 
defined by 



P* = /3*(d,£) 



if^ = (Regime 0), 

if^(0,d,£)>0 (Regime I), 

Pii if 1^(0, d, e) < and (^5(7*, d, e) > (Regime II), 

7* if ^p{-i*,d, e) < (Regime III). 



(4.25) 



In this definition, /?// = j3ii{d,e) is the unique solution of ip{-,d,e) = with ip defined by 

ip{P,d,e) := -det(A£*(/3,d,e)) = -dct[a(/3, d)-ie(e)j , (4.26) 
Ae* = A£*(/3,d,£) :=£*(/?, d,£)-eJ(/3, d,e), 
and the yet undefined functions arc specified below. 

The four regimes have the following crystallographic interpretation, which follows from the 
construction of the optimal microstructure in the calculation of W. 

Regime 0: The material is homogeneous and the energy does not depend on the microstruc- 
ture. This occurs when a2(e^ ~ s) ~ Oii{eJ — e) = 0. 
Regime I: There exist two optimal rank-I laminates. 
This is characterized by 

det [{d2ai + ^102)"^ (a2(e^ - e) ~ ai{ej - e))] < 0. 

Regime II: The unique optimal microstructure is a rank-I laminate. 
This regime occurs when the function 

[0, 7*] 9 /? ^ det [(daai + dia2 - l3Ty' ("2(4 " e) " "ilef ^ e)) 
has a unique root (which we denote by /3//). 

Regime III: There exist two optimal rank-II laminates. This regime is present if the 
operator (d2ai + dia2 — "f*T) is invertible and 



det 



(d2ai -f- dia2 ~ j*T) ^ (a2(£^ - e) - aiCef ~ e)) 



> 0. 



For illustration, we visualize prototypes of rank-I and rank-II laminates in Figures [3] and S) 



To complete the definition (|4.26p . set 

a(/3*,d) 
eie) := 
e*=e*{l3*,d,e) := 
ei(/3*,d,£) := 
e2(r,rf,£) 



d2ai + dia2 - (3*T, 

a2{el - e) - ai(ef - e), 

a'\(i\d)e,{P\d,e), 

(a2 - P*T)£ - d2{a2e^ - aie{), 

(ai - /3*T)£ + di{a2e^ ~ aief ). 



(4.27) 
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Figure 3: A two-phase rank-I laminate in two space dimensions with corresponding normal 
vector. The strains arc constant in the shaded and in the unshaded regions. The volume 
fraction of both phases, 0.5 in the picture, is prescribed by the macroscopic parameter d. 




hi h2 



Figure 4: A two-phase rank-II laminate in two space dimensions. The widths hi and /12 of 
the slabs should be much larger than the thickness of the layers between the slab. 



Kencc 62 — el ~ [a{/3* , d)] ^e{e). 

Explicit computations of ddW and d^W arc lengthy. We recall the following results for the 
partial derivatives W 'm D = 2 which are proved in [3 . Set 

cr^ := ade, - e, ), 
a := diCTi + d2(J2 ■ 

Lemma 1. Let D ~ 2 and let ai and T commute. Then 
dW , 



^^(d,e) = diai {el - ej) + ^202 (^2 - ^2) 

^ I 7*c?ic?2Q!~^(7*, d){ai — a2)T (ej — e*) in Regime III 
I else. 



— (d, e) = di (a2 - 7*T)a-i (7* , d)ai {el (7* , d, e) - ef) 



(4.28) 



Alternatively, in Regime III, 
dW , 

d,e) = di[a2 - Y'J-')oi "(7 

-I- d2[ai - YT)a'\^*,d)a2 (£^(7*, e) - el) . (4.29) 
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Moreover it holds 



^(d,e) = a* -.As* + Wi{el) - W2{e;) + <j /3*did2^|| A£*f in Regime II, 



in Regimes and I, 

P*did2^\\Ae*f in Regime II, 
{di — d2)j*(p{Ae*) in Regime III. 



(4.30) 



^ = J ((d2ai+diQ2-/3T)-i(TAe*)):TAe* -TtCf/imc 11, (4.31) 

1 otherwise. 



(4.32) 



^ j iid,c.,+d,c.,-pT]-^iTAe')y-TAe' ("i " "2)a 'TAe* m Regime II, 
de 1 otherwise. 

4.2 The scalar setting for D = 3 

The scalar setting in three dimensions is characterized by the ansatz 

Xl 

u{xi,X2,X3) := I X2 ] (4.33) 

r]{xi,X2) 

for the deformation, where 77 is a scalar function (hence the name 'scalar' theory, see [S], 
[5]). Physically, (|4.33p corresponds to anti-plane shear in the xs-plane. Since 

e(u) = 1 dx2v\ , 
\dxiri dx2V / 

V?7 determines the strain. This justifies to work with vectors / = V77 G M^, not with 
matrices e(u). We replace p.l3p by 

W{d,f):= inf inf -f diWi{S/fi) + d2W2{S/fi) dx, dG[0,l], (4.34) 

{d)=d V\an'=f-x Jn' 

Je{o,i} 

where the domain of integration Vl' is now two-dimensional. Therefore, d £ BV{il.'; {0, 1}), 
and it is well-understood from the context that 



denotes the two-dimensional average. In addition, in (|4.34p we adapted the common notation 
for the scalar case and wrote / instead of e for the strain. Finally, we set 

wdf) ■■=\o^.{f'fI)-if-fT) + w, 
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with ff the transformation strains in the scalar setting and positive definite matrices. 

In [T3|, partial results for W in the non-scalar, i.e. vectorial setting are available for D = Z. 
Yet, in general, the computation of ddW , d^W for I? = 3 as required for assumption {A) 
remains currently open due to its complexity. This is the reason why we restrict ourselves to 
a special three-dimensional case. Our main result in this subsection is the following theorem 
that provides us with an explicit formula for W in the scalar case. 

Theorem 3 (Representation formula for W and D = 3 in the scalar setting). The functional 
W , given by ^.34^ , satisfies the explicit representation formula 

W{dJ) = diWiiVT]l)+d2W2iV7j;) (4.35) 

with 

V77i*(d, /) (dsai + diaa)"' [aa/ - ^2(02/2^ - ai/f )] , (4.36) 

V7?2*(d,/) ■■= id2ai + dia2)~'[aif + diia^f^ - a^fl)]. (4.37) 



Theorem |3] states in particular that the scalar elastic theory in 3D is directly related to the 
non-scalar geometrically linear elasticity theory in 2D which was discussed in Subsection l4.1l 

Proof: We apply the translation method, see for instance (TB], [H] for an overview. In 
general, the function ip which describes the translation is only quasi-convex. Here, as a 
benefit of the scalar theory, we can pick 1^9 as a convex function, cf. (|4.4ip . 
Starting from the inequality 



'/'(/) ^ inf + (p(Vt]) dx 



which is satisfied by any quasi-convex function (p, we obtain, similar to the reasoning in 
for the non-scalar case in D = 2 dimensions. 



W{d, f) > max min inf 

I3>0 deBV{n':{0,l}) v\dn'=f-x 

ly,-/?.^ convex ^^^^^ 



(4.38) 



Consider 



Then 



^-..di V?7 dx ~ -f?^,d2Vri dx 

Vm{d,v):= ^'': ~ ' , Vmid,v)--- Z J . • (4.39) 
j-^,di dx fj^,d2 da; 



diVrjiil ri) + d2\/r,2id, ri) = /. (4.40) 



Next we apply Jensen's inequality which is possible since Wi — /Stp is convex for i = 1,2. 
Furthermore we possibly enlarge the set of admissible functions in the minimization over rj 
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in (|4.38p by considering the set of admissible functions V771 , V?72 G IR^ with the constraint 
glOl). Then (|4:38l) becomes 

WidJ)> max min \ diWi{Vm) + d2W2{S/7^2) + Mf)} ■ 

VKi-^.^ convex (iiVl)l+rf2Vl72=/ 

If I? = 2, (p is chosen to be — det(e(u)). Here we set analogously 

fiWr,) (9,77)' + {dyvY > 0. (4.41) 

Since ip is quadratic, there exists a unique linear operator T : R'^ — )- such that (p{f) = 
\Tf:f, here simply Tf = 2/. It holds 

2 

i=l 

2 

= ^ ( V?7, ) + /3 [ y(diV7?i + d2Vm) - diLp{Vrii)-d2Lp{Vm\ ] ■ 

By a direct computation the remainder term R can be rewritten, 

2R = T{diWiu + d2 Vrya) : (di V?7i + ^2 V772) - rfiTV?7i : V771 - d2TVT]2 : 

= di<i2 [ - TWm : V?7i - TV?72 : Vr72 + TV?/i : V?/2 + TWm : V?7i] 

= -did2T(V772 -V77i):(V?72- V771) 

= -2did2</5(V?72 - V771). 

So we have found 

W{dJ)> max min |diiyi(V77i) + d2l1^2(V772) - /3rfirf2'p(V?72- Vt^i)}. 

/3>0 Vr)i.Vi)2eR^ J 

VVi-ZJ.^ convex diVrii+d2Vri2=f 

(4.42) 

Next we compute the optimal strains Vr^jf, Vr^j- After differentiating the argument on the 
right in (|4.42p for fixed /3, we obtain 

ai(V77i* - fl) - «2(V?;2* - fl) + /^^(Vr;* - Vr;*) = 0. (4.43) 

Using the constraint diVijl + d2V772 = /, after rearrangement this yields the formulas 

(dia2 + d2ai - l3T)V'ql = {a2 - (3T)f - ^2(^2/^ - aifi)-, 
(dia2 + d2ai - /3T)V77* = («i - /3T)/ + ^1(0-2/2^ - a^fl). 

Setting /? = (see below), this coincides with (|4.36p . (|4.37p . 

The maximum over /3 in (|4.42p is attained at /3* = as did2^{'^ri2 — '^rji) > 0. Hence, the 
translation is trivial in this setting. We obtain 

W{dJ) > diM^i(V?7i*) + d2l^2(V?72*) =: W-{dJ). (4.44) 
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It remains to estimate W from above which amounts to showing that V77J', V772 yield an 
optimal microstructure. Plugging in any microstructure d with (d) = d, from the definition 
(|4.34p of W, we get the upper bound 

W{dJ) < inf / diWi{\7T]) + d2W2{\7v) <ix W+{dJ). 

As the domain of integration ft' in the scalar case with D = 3is two-dimensional, finding the 
optimal microstructure subsumes to the non-scalar case in D ~ 2 dimensions. Hence, as can 
be seen from (|4.25|) . /3* = occurs either for Regime (no microstructure) or Regime I. The 
latter is equivalent to (p{e2 ~ £*) > 0. From a well-known argument, see [Ml Lemma 4.1], 
this implies that ej, 62 are compatible. Besides this compatibility condition, any optimal 
microstructure must also satisfy the equilibrium condition 

[[a]]n = (4.45) 

which states that the jump of the stress in the normal direction of the laminate must vanish. 
Equation ()4.45p originates from the Euler-Lagrange equation of the variational problem. 
Here, it is automatically satisfied since, due to (|4.43p . [[(j]]n = crj — o-J* = 0. 
So there exists a unique rank-I lamination microstructure (unique up to a sign ±1 of direction 
of the laminate) that connects e* and 62, and the strain of phase i is Vi]*, i = 1,2. This 
microstructure is optimal, W+{d, /) = W-{d, /), which proves (|4.35p . □ 



4.3 The case D = 1 



The methods of the previous subsection also apply to derive rigorously an explicit expression 
for in D = 1. Here, quasi-convex functions are always convex, and the vector and scalar 
settings coincide. However, we use the notation in analogy to the vector case. 

Theorem 4 (Representation formula for in = 1). For D — \, the functional W , given 
by h2.13\) . satisfies the explicit representation formula 



with 



W{d,e) = dWi{el) + {I - d)W2{e*2) 
02 (e — d2e2) -\- d2aie\ 



e\{d,£) 
e;(di,£) 



d2ai + dia2 
ai{e — diej) + dia2£^ 



d20Li + dia2 



The partial derivatives of W are given by 

d2{ 



si) 



4) 



dW , , , aia2[di{e 

^ ' 

oe d2ai + aiQ!2 

^(d,e) ^ W,{sl) W2{e*2) + 
od [d2ai + dia2) 



(4.46) 

(4.47) 
(4.48) 

(4.49) 



(ai - Q;2)e^ + diaiie'I)'^ - d2a2(eD^ 
+ {a2e2-aiel + (a2 - ai ) (dief -|- d2e^)) £ 



{d2ai - dia2)e-^ e 



(4.50) 
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Proof: We apply the same methods as in the proof of Theorem [3] 

First, we estimate Wi, W2 from below by their convexification. Similar to (|4.42p . we have 
for any convex function ip and /? > 0, 

W{d,e)> max min {d{Wi - Pip){ei) + {I - d){W2 - M(^2) + M^)] ■ 

M/.-/3<^ convex dei + (l-(i)e2=e 



Secondly, we choose the optimal /3 and (p. The later is an art when D > 1. With arguments 
identical to the scalar case in D = 3, it holds (3 = 0, leaving us with the estimate 

Wid,e)> min {dM^i(ei) + (1 - d)M^2(£2)} 

dei + {l — d)s2—€ 

Next we compute the optimal strains Si, £2 on the right. For = 1, only £1 = e is 
admissible and hence optimal. Now, let < < 1. We can then resolve the constraint by 
setting £2 = ^j-rj' and need to calculate 

imn|dM^i(£i) + (l-d)M/2(^^^)|. (4.51) 

After differentiation with respect to £1 and using (|2.12l) for the optimal £* in (|4.51[) . we are 
left with 

del 



ai(£t -£i ) = Q2(^ ^ _^ - £2 

Rearrangement and resolution of del + (1 ^ d,)e2 = £ gives the formulas (|4.47p . (|4.48p of 
the optimal strains. Since £*(1,£) = £, these formulas also state the correct solution when 
d=l. _ 
In the final step, due to the definition of W, it holds 

W{d,e)< inf -f dWi{e{uj) + {1 ~ d)W2{e{u)) dx =: W+{d,e) 

and d e BV{il] {0, 1}) may represent any microstructure with (d) = d. The explicit con- 
struction of the optimal microstructure for which the upper bound and the lower bound 
coincide, 

W-id,e) = W{d,e) = W+id,e), (4.52) 

is again as in the proof of Theorem [3] This leads to (j4.46p . 
Finally, the verification of (|4.49p , (|4.50p follows from 

'^^^^^^^1 — r:^ — ' — n — ^^-^^^ 

oe d2ai + aia2 ue d2ai + dia2 

and 

del o;2 [a2£2' — '^1^1 + ("1 ^ a2)£] 



dd (^2^1 + dia2)^ 

de2 _ Oil [a2e2 - aiel + (ai - Q!2)£] 

dd {d2ai + dia2Y 



(4.54) 
(4.55) 
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which all can be verified by elementary computations. 

Using the relationship c?ie* + ^2^2 = e we find (|4.49p and finally 

8W rlr* 8r* 

— = W,{el) - W2[el) + diW[{el)-^ + dMe*^)^, (4.56) 

= Wi{el)-W2[el) 
^ aia2{a2£l ~ aief + (qi - a2)e) [dije'l - el) + ^2(62 - eT^)] 

(d2ai + dia2Y 

Using 

- ) + ^2(^2 - ef) = £ - c'2£i' - rfief, (4.57) 
this simplifies to identity ()4.50p . □ 

5 Existence results for the AC-CH system with mi- 
crostructural energy densities 

Above we have already remarked that the energy Wiin satisfies assumption {A) and thus 
Theorem [1] applies. In this section we investigate whether the results transfer to the geo- 
metrically linear theory of elasticity stated in the previous section. 
The following statement is a consequence of Theorem [1] 

Theorem 5 (Existence of solutions for geometrically linear elastic energy). Let W ^ W 

he a function on [0,1] x R^J^^. Assume that the mobility tensor M he positive definite. 

Moreover, let ip he given hy ^ and the initial data (ao,bo) satisfy ilO\}. 

In D = 3, let W he given hy ^.35^ corresponding to the scalar setting. If D = 2, let 

(CI) P*{d,e) he independent of e, and 

(C2) Ui and T commute whenever 13* S {7*,/?//}. 

If D — 1, let W he given by {4. 40) . Then there exists a solution (a, 6, u) to (I2P^i3i that 
satisfies 

(i) a,beC°^h{[Q,J]-L^{n)), 

(11) dta, dtbeL^inr), 

(lii) ueL°° (0,r; H^{n; M^)), 

(iv) The feasible parameter range of {a,b) is given hy hlO\}. 

Proof: We first consider the case < 3. The idea is to apply Theorem [1] Therefore we 

1. verify that W satisfies assumption [A) for d G [0, 1], 

2. extend to R x Mf^^,^ in such a way that {A) is fulfilled for all d e M. 

Then Theorem[T]can be applied and we obtain by part (iv) therein that effectively d € [0, 1] 
in the system ([U-([3]). Hence Theorem [5] follows with yV = W &s asserted. 
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We start with the second step and give the proofs of the former in dimensions D = 1, 
D = 2, and the scalar case in D = 3 thereafter. We introduce the foUowing extension 
Wo(-,e) :M^M of M?(-,e) by 



' W{d,e), 0<d<l, 

gi{d,e), 0> d> -1, 

-d+l + W{0,e), d<-l, (5.58) 

&(d,e), l<d<2, 

I rf- l + t?(l,e), d>2 



for functions gi, g2 that we may choose such that Wa{-, •) G C^(M x R^j^^^) and that Wq 
fulfils {Al)-iA2) for d < and d > 1. 

We further mention that the function Wq defined in (|5.58p is positive and coercive for 
\d\ — >■ oo. We note that this construction is required for the first part of the proof of 
Theorem [1] where polynomial approximations of -0 are constructed, where d ^ [0, 1] may 
occur. 

It remains to show that W fulfils (A) for d e [0, 1], which we prove independently for D = 1, 
D = 2, and the scalar case in D = 3. 
(i) Verification of (A) for D = 1: 

From (Pl^ . glU), (|071) and gUl), h follows W e C^{m. x R). 
From (|449)) . we obtain 

2 

2—1 

The condition (^1), restricted to d e [0, 1], follows, since by (jimi . (f05|) for z = 1, 2, 

"^^^^^ > by assumption. From (lOT)) . for i = 1, 2 and d e [0, 1] wc have 

by similar arguments 

\eUd,e)\<c{\d\ + \s\ + l). (5.59) 

This leads to {A2)i, since by 

|W^(d,£)| < 11^1 (£l)| + |M/2(e;)| < c (|etp + + l) 
<c(Mp + |£p + l). 

For d G [0, 1] we easily compute, for j = 1, 2, 

\d,e*id,e)\<ci\e\ + l). 

We use this estimate in ()4.56p . to obtain by (|5.59p 

\ddW{d,s)\ <c{\sUd,e)\^ + \e*M,e)\' + 1 + (|£^(d,e)| + \e;id,em\s\ + 1)) 
<c(Mp + |ep + l) 
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which is (v42)2 restricted to d S [0, 1]. 

For d e [0, 1] we derive with the help of gH]), (|437l) . ((539| 



\deWid,e)\ < 



{\e~e{\ 



min{ai ,a2 } 

< max{ai,a2} c{\e\ + 1) 
<c(|£| + l), 

showing the vahdity of (-42)3 restricted to c? e [0, 1]. 
(ii) Verification of (A) for D = 2: 

The regularity of W required for (A) follows directly from (|4.23p and the related definitions. 
The constant 7* solely depends on ai, a2 since 7* = min{7i,72} and ji is the reciprocal of 
the largest eigenvalue of a^^^^Ta^^/^ see g^ll). Since e [0,7*], 



W*id,e)\<c 

for a constant c independent of d and e. 
From ([ETO)) . tfn7\ we obtain for d S [0, 1] 

|£*(/3*,d,e)|<c(|d| + |£| + l). 

With (12121) and ((^^ this shows for d € [0, 1] 



(5.60) 



(5.61) 



|/3*did2det(£*(d,£) - el{d,e))\ < c{\el{d,e)\^ + |e^(d,£)P + l) 



<c(|dP + |£p + l) 



From this we easily verify (-42) 1 restricted to c? e [0, 1], as the first two terms in (j4.23p can 

be estimated as in the one-dimensional case. 

In order to show (-42)2, we apply (|4.30p and (|5.59p to find 



r(d,£)i 



^d,a,;(£* -£f) 
1=1 

1), 



<c(|d| + |£| 
|A£*(ci,£)| = |£;-£j| 

<c(|d| + |£| + l), 

|^(A£* (d, £))| = I det(A£*(d, £))| < |£*(d, £) - elid, £)|2 

<c(|dp + |£p + l). 

Finally, from (|4.3ip and the fact that T is an isomctry w.r.t. the Frobcnius norm, 

d/3* 



dd 



id,s) 



< c. 



The terms Wi{e*{d,e)), i = 1,2 remaining in (|4.30p can be estimated as in the one- 
dimensional case. This verifies (-42)2 restricted to d e [0, 1]. 
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Now we prove (.42)3. For the regimes 0, I and II this foUows directly from (|5.59p . When in 
Regime III, we start from (|4?29| and find with ([53T|) for d e [0,1] 

\di{a2 - i*T)a-^a^{el - )| < c{\e\\ + 1) < c {\d\ + \e\ + 1) . 

The vahdation of (.41) relies on (|4.28p . By assumption (i) of the theorem, = and, as 
d £ [0, 1] is fixed in (.41), is a constant tensor. Assmnption (.41) then follows from 

(£^(^,£2) - el{d,ei)) : (£2 - £1) = a-\a2 - p*T){e2 ~ £1) : (£2 - £1), 

a similar equality for (£2(rf, £2) ~ £2(^7 ^i)) • (£2 ~ £1), and the positive dcfiniteness of 

and (a; — P*T), i — 1,2. Indeed, the latter is equivalent to Id — l3*a.^ ^^"^Ta.^ positive 

1/2 1/2 

definite, and since /3* < 7*, this is implied by all eigenvalues of Id — 7*0, Ta^ being 
positive. However, from (|4.24p and the definition of 7^, this is true. 

For the scalar case in 13 = 3, comparing ()4.23p and (j4.35p and since the partial derivatives 
coincide, the theory subsumes to the geometrically linear theory for D ^ 2 and /?* ~ 
0. Hence we can proceed as above and show that W given by (|4.35|) satisfies (.4). The 
assumptions (CI), (C2) are not needed here because /?* = 0. □ 

The following theorem states the uniqueness of weak solutions for the system with geomet- 
rically linear theory. 

Theorem 6 (Uniqueness of weak solutions). Let W = W be given by {4.35\j for D = 3 
(scalar setting), by ^.23\ l for D ~ 2, or by ^J^J^ for D = 1. Assume further that M = 1 and 
the elastic moduli of the two phases are equal, i.e. oi\ = a^. For D = 2, let det{£ 2 —ej) < 0. 
Then the solution (a,5, u) in Theorem\^is unique in the spaces stated there. 

Proof: (i) The case D = 1. 

For ai =0:2, the equations (|4.47|) . (|4.48p read 

el{d,e) = £ + d2{ei - £^), £2(rf, £) = £ + di[e2 - e{) 
which implies £2 ~ £2" = £* ~ £1" leading to Wi(£*) — ifi — M^2(£2) ^ ^^'2- This yields 

W{d, e) = ^(£^ - elY + diwi + d2W2 = + di{e2 - £f ) - £2")^ + diWi + d2W2. 
In this special case, the partial derivatives of W are 

ai(e + di{s2 - e'{) - £^)(£j - £f ) +wi -W2, 

a^{e + d,is^ ~ej)~sl). 

Now we pass through the steps in the proof of Theorem [2] and make the necessary modi- 
fications. Let again {a^ ,b^ ,u^), fc = 1,2 be two solutions to (HI)-© and set a := — a^, 
b:=b^ - b^, d^ := + b^\e^ := £(u'=) for k = 1,2 anA d := d^ - d^ = a + b, e -.^ ~ . 



dW 



{d,e) 
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The equations ((XTO)) . (j??^ now read 

J [dtbr] + AV6 • Vry + ai (e + d(4 " )) (e^ " ^Dv] dx dt 

^ (g'ia^ + b^) - g'ia' + b') + g'{a'-b') - g'{a^-b^)) v + ^K^brj 



Or 



Or 



J ai{e + d{e2 -ef))e{u)dxdt^O. 



dxdt, (5.62) 
(5.63) 



Analogous to ()3.22p we find 

^Hto)\\M+ I A|Vap+«i(e + d(e^-ef))a(e^-£f))dxdt < 



2Kia 



Choosing 77 
obtain 

A 



g'{a^+b^)- g'{a'^ + b^)\ + \g'{a^-b^)- g'{a^-b^)\ \a\dxdt. (5.64) 
^ — &^)A:'(o,to) in ()5.62p and adding the resulting equation to ()5.64p . we 

dxdt (5.65) 



2ll«(^o)|pM + ^||fo(MllL^+ / A(|Vap + |V6p)+ai(£ + d(£^-£f))' 



< / 2(Ki|a|^ + «:2|6|')+ / ^ [|.g'(a'+*')-ff'(a'+6')| + |5'(a -&')-ff'(a'-^'')|] + l&l) di. 

It is noteworthy that we are only able to estimate j^^^ ai(e + (i(e|^ — e|"))^ dx dt whieh differs 
from the total mechanical energy /^^^ W{d, e) dx dt. 

With the Lipschitz continuity of p', and by applying the inequalities of Young and Gronwall 
to (|5.65p , we obtain as in the proof of Theorem [5] 

ai{e{u)y dxdt = 0. 
Or 

With Korn's inequality, m = in ^7-. 
(ii) The scalar case in D = 3. 
By Theorem [3j for ai = 02, 

V7?2*(d,/) =/ + dl(/2^-/f). 

Consequently, Vr/J' - = Vr/2 - /J leading to WiCV-ql) - wi = VK2(V77|) - W2. Since 
/3* = 0, the formula for W has now the same structure as for D = 1 and the proof follows 
exactly as in the first part (i). 
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(iii) The case D = 2. 

By assumption, ip{/3,d,e) > 0, which resuhs in /3* = in O7-. In addition, for ai ~ a2, by 
Eqn. g^lZl), 

slid, e) = a^' {{a, - /3*T)e - d2«i(4 - )) , 
£;(d, e) = ar' (("1 - l3*T)e + d,a,iel - e{)) . 

Consequently, 

d,e) = - det((e; - £l)(d, e)) = - dct(4 - e^). 

So again el — ef = £2 ~ ^2 ^ ^1(^1 ) — wi = W2(£2) ~ "^2, and the proof can be carried out 
as before. □ 



6 Conclusion 

In this article we derived extensions of the Allen-Cahn/Cahn-Hilliard system to elastic 
materials. In particular we included (i) the linear elastic energy derived by Eshelby, (ii) a 
geometrically linear theory of elasticity for Z) < 2, and (iii) a geometrically linear theory of 
elasticity for Z) = 3 in the scalar setting. For all three cases we showed in a mathematically 
rigorous way the existence and the uniqueness of weak solutions, asserting the correctness 
of our approach. As a future goal, it is desirable to generalize the proof to the non-scalar, 
i.e. vectorial case in D = 3 and to Z? > 3. 

The generalized AC-CH models contain as special cases both the Allen-Cahn and the 
Calm-Hilliard equation [211 with linear elasticity. In our work the existence and uniqueness 
of weak solutions to the AC-CH model with linear elasticity follows in a straightforward 
way since Assumption (A) required in the existence result Theorem [1] can easily be checked 
to hold true for the linear energy. This had not been shown before. 

We point out that for the new cases (ii) and (iii), the formulas collected in Section [3] for the 
specific energy W are essential for any numerical studies of the AC-CH models extended 
to microstructure. For related investigations and numerical methods we refer the interested 
reader to [7| and [8]. 

Besides the limiting assumption of constant temperature (cf. Introduction), the most im- 
portant pending restriction is the postulation of small strain, included in That is, it 
would be desirable to combine the AC-CH system with energies within the framework of the 
geometrically non-linear theory of elasticity. The technical problems of a large strain theory 
are striking, see, e.g., [H [5l I15[ 125] . and it is in many cases not known how to compute 
explicit formulas for the relaxed energy functionals. 

Similar mathematical difficulties arise when one wants to extend the AC-CH systems to 
elastic materials which are of polycrystalline structure. This is an open problem for future 
research. For related modelling aspects we refer to [TU], for related mathematical aspects to 
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